Power spectrum sensitivity of raster-scanned CMB experiments 
in the presence of 1// noise 

Tom Crawforcj3 

Department of Astronomy & Astrophysics and Kavli Institute for 
Cosmological Physics, University of Chicago, Chicago, IL 60637 

We investigate the effects of f// noise on the abihty of a particular class of Cosmic Microwave 
Background experiments to measure the angular power spectrum of temperature anisotropy. We 
concentrate on experiments that operate primarily in raster-scan mode and develop formalism that 
allows us to calculate analytically the effect of 1// noise on power spectrum sensitivity for this class 
of experiments and determine the benefits of raster-scanning at different angles relative to the sky 
field versus scanning at only a single angle (cross-linking versus not cross- linking) . We find that the 
sensitivity of such experiments in the presence of 1// noise is not significantly degraded at moderate 
spatial scales {£ ~ 100) for reasonable values of scan speed and 1/f knee. We further find that the 
difference between cross-linked and non-cross-linked experiments is small in all cases and that the 
non-cross- linked experiments are preferred from a raw sensitivity standpoint in the noise-dominated 
regime — i.e., in experiments in which the instrument noise is greater than the sample variance 
of the target power spectrum at the scales of interest. This analysis does not take into account 
systematic effects. 



I. INTRODUCTION 

Scanning strategy is an important component in 
the planning of Cosmic Microwave Background (CMB) 
anisotropy measurements, particularly in the presence of 
long-timescale drifts in the gain or DC level of the in- 
strument response to incoming signal. These drifts may 
be due to intrinsic quantum processes in the detector or 
readout electronics, temperature changes in the instru- 
ment environment, or slow changes in the behavior of an 
external noise source, to name but a few examples. It 
is common to lump all of these under the heading 1/f 
noise, because the spectral density of such drifts often 
approximates a power law in frequency. The effect of 
such drifts on instrument sensitivity to sky signals of in- 
terest depends on — and often drives the design of — 
the instrument scan strategy. 

It has become conventional wisdom in the field that 
sensitivity in the presence of 1// noise depends crucially 
on the degree of cross-linking in the scan strategy — that 
is, the number of different scan angles from which each 
pixel in the map is observed [e.g., [ll, 0]- Cross-linking is 
also often invoked as a means to minimize the effects 
of systematic effects such as scan-synchronous ground 
pickup and chopping mirror offsets; in this work, how- 
ever, we only address cross-linking as a means to combat 
the effects of 1// noise. By the criterion of degree of 
cross-linking, the raster-scan strategy — in which the in- 
strument beam is scanned back and forth across the sky 
field in azimuth and stepped in elevation — is maximally 
sub-optimal, in that the majority of map pixels are only 
ever observed from one scan angle. However, for many 
sub-orbital (i.e., ground- and balloon-based) CMB mis- 
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sions, this is the preferred mode of observing, because 
any elevation component to a scan will result in a large 
scan-synchronous signal from the changing optical depth 
of the atmosphere. For most sub-orbital platforms, some 
degree of cross- linking can be achieved with azimuth-only 
scans by observing a sky field at different times of day, 
using sky rotation to change the scan angle. For obvious 
reasons, this technique does not work for instruments 
observing from the South Pole. Despite this, successful 
measurements of CMB anisotropy have been made from 
the South Pole using single-dish instruments in raster- 
scan mode , and instruments at the Pole that comprise 
a major part of the present and near future of CMB sci- 
ence are currently operating or plan to operate primarily 
in raster-scan mode. [1, H, Q 

It seems therefore worthwhile to investigate the lim- 
its on CMB power spectrum sensitivity from a raster- 
scanned instrument in the presence of 1/f noise. While 
the effects of 1// noise on CMB power-spectrum sensi- 
tivity have been investigated by other authors, (e.g., Q), 
the specific case of raster-scanning instruments has not 
been dealt with since the treatment of [2] , the conclusions 
of which are the focus of this work. In particular, we wish 
to determine the relative merits of a scan strategy with 
a single scan angle versus one in which the scan angle 
is changed on the timescale of the rotation of the Earth, 
as these are the maximum amounts of cross-linking one 
can achieve with azimuth-only scans from a Polar and 
mid-latitude platform, respectively. 

To achieve this, we first develop formalism that en- 
ables us to calculate analytically the effect of 1// noise 
on power spectrum sensitivity for raster-scanned exper- 
iments. This effort is similar to that of 8], who per- 
formed the calculation for instruments that scan along 
interleaved rings. We apply the formalism we develop 
to two fiducial scan strategies, one with no cross-linking 
and one with the minimal cross-linking available from a 
sub-orbital platform in the mid-latitudes, using two pa- 
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rameterizations of 1// noise: a toy model which provides 
a simple instructive result and a more realistic model. Fi- 
nally, we discuss the results of these calculations in the 
context of earlier work on the subject, particularly that 
of i. 



II. FORMALISM AND TOY-MODEL RESULTS 

To try and get an analytical feel for the effects of 1// 
noise on CMB power spectrum estimation from raster- 
scanned observations of a square sky field using a single- 
pixel instrument (the generalization to an array instru- 
ment is trivial if the 1// noise is uncorrelated between 
array elements), we will begin by modeling 1// noise as 
a step function with one white-noise value above a par- 
ticular frequency and another, much larger white-noise 
value below this frequency. After extracting the simple 
but illustrative result from this case, we will apply the 
calculation to a slightly more realistic parameterization 
of 1// noise. 

A. Sensitivity from a single observation 

First, we note that for a small enough sky field, (<20° 
in each direction), we can use the flat-sky approxima- 
tion to decompose the CMB fluctuations on this field in 
Fourier modes. Following 

AT{x, y) = (27r)~2 j (fk a(fc)e'('^-^+'=''^). (1) 

The a{k) are zero-mean Gaussian variables with variance 

{a{k) a*{k')) = P{k)5{k-k'), (2) 

where P{k) ~ Ci and A: — |fe| ~ ^ in this approximation. 

We wish to estimate P{k) from our (noisy) map. Fol- 
lowing we note that for a set of N independently 
measured modes with |fe| = k and variance P{k) + 
pnoiscj^^-j^ the variance on our estimate of P{k) is: 

[dPik)]" = {{P^^\k)-P{k)f) (3) 

= |(P(fc)+P--^(fc))'. 

Our sensitivity to the amplitude of an arbitrary spatial 
mode on the sky Ar(xj,, yp) = a{mp) nip is given by: 

aMaK)) ^^°^<^'-^'"^^'' 

where Wppi is the pixel-pixel weight matrix of our map.^ 
If we have made a minimum- variance map from our ob- 
servations, then 



Wpp, = AjpWtt'Afp', 



(5) 



^ This is only strictly true for modes that are perfectly resolved by 
the instrument beam. 



where Atp is the pointing matrix (the operator that de- 
projects a spatial mode into the timestream), and Wu' 
is the inverse of the timestream noise correlation matrix 



Wu 



{n{t}n^{t')y 



(6) 



If we have 1/f noise in our timestream, then Wu' will 
not be diagonal, but if the properties of the noise are 
stationary in time, then Wu' is circulant, and its Fourier 
transform 



Wff = Ff.Wu'Fj,,, 



(7) 
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is diagonal. We can then rewrite equation |4] using equa- 
tion [5] and inserting two pairs of forward and inverse 
Fourier transforms: 



Wa 



nip {AfpWu'At'p') nip' 

{Atpnipf F^FWu'F^F [At^p^nip,) 

^WFAm] (/)p 



(8) 



/ ^noisc 

/ 



(/) 



where [FAm\ (/) is the Fourier coefficient at temporal 
frequency / of the time-domain function that results from 
deprojecting mode ni into the timestream. (In other 
words, we have derived the unsurprising result that if 
the timestream noise is uncorrelated in the Fourier do- 
main, our sensitivity to a mode on the sky is the sum of 
the (squared) Fourier coefficients of its deprojection into 
the timestream, weighted by the inverse Fourier-domain 
variance.) 

For a raster-scan observation beginning at the origin of 
our sky field, scanning at angular speed Vs at an angle 9^ 
from the x direction of our field, and stepping a distance 
Ay in the y direction of our field every seconds, a 
given Fourier mode from the expansion in equation [1] will 
deproject into the timestream as: 



AtpATp{k) = a{k) X 



(9) 



^i[kxVs cosBs t~\-ky{vs sin 0g t+Ay floor(f/Tr))] 



i[kx cos 8s + ky {vs sin f?g+Ay/ri-)]i 



a{k)e 



where the "fioor" function returns the nearest integer 
less than or equal to its argument. The second line of 
equation[5]approximates the stepping in ?/ as a continuous 
(very slow) scan in y. It is clear that the deprojection 
of a single Fourier mode into the timestream of a raster- 
scanning experiment has power at only a single temporal 
frequency: 



27r 



— 1 Vs sm f 
Zn 



Ay 



(10) 



In the absence of mode-mode correlations, the vari- 
ance on our measurement of the coefficient a{k) is sim- 
ply proportional to the timestream noise variance at the 
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frequency /(fc): 

(a(fc)) 



Wa{k) 



(11) 



— P"-- (/(fe)) . 



where Ns is the number of times the instrument is 
scanned across the field between elevation steps. As tra- 
ditionally defined [e.g. 0], noise with a 1// component 
has a power spectrum: 



= (/)=< 1 + 



/i 



knee 



(12) 



where tr^ is the variance in the high-frequency "white" 
part of the power spectrum, and /knee is the frequency 
at which pnoise^y^ 2(7^. If we approximate this as 



pn 



_ al, if / > /kn 



oo, otherwise; 
then the variance on our measurement of a{k) is: 



a' (a(fc)) = 



_fcTl/N,, if/(fc)>/k„ 



oo, 



otherwise, 



(13) 



(14) 



and the variance on our estimate of the CMB power spec- 
trum P{k) from this map will be 



[SPik)]'- 



Ngood{k) 



{P{k)+al/Nsy 



(15) 



where iVgood(fc) is the number of modes with |A;| = k 
that satisfy f{k) > /knee. 

If we define 9k as the angle between the x direction of 
our sky field and the direction along which a particular 
Fourier mode is oscillating, then 



= |fc|cos6'fc 
ky = |fc|sin6'fe. 



(16) 



and our criterion for goodness can be expressed as: 

^Tr/knee 



COS ( ffc 



smek — - 

VsTr 



> 



(17) 



For raster timescales Tj- that are sufficiently long com- 
pared to l//knco, we can neglect the raster term in equa- 
tion [T7] and write: 



cos Wk 



s)\> 



^Tryknee 

kv. 



This allows us to define a critical angle 

e^ik) = cos-i (^^^) , < 0c < 



(18) 



(19) 



If 27r/knee/(fc^^s) > 1, then there are no good modes at 
that k value, and Oc{k) = 0. (Note that for observations 



at a single scan angle, the choice of coordinate system 
that defines k^ and ky is arbitrary, in which case we can 
define them such that 6*5 = 0° , and our criterion for good- 
ness is then simply k^ > 27r/knee/i's-) 

The fraction of good modes with a given |fc| = fc in the 
sky map can now be expressed as: 



-/Vtot(fc) TT 



(20) 



A^tot, the number of independent modes within a bin of 
width Afc centered on k for a map covering a fraction 
/sky of the sky is given by: 



iVtot(fc) = 2fc/skyAfc, 



(21) 



which allows us to write the variance on our measurement 
of P{k) from this observation of this patch of sky as: 



[5P{k)\ 



2fc/skyAA: 20e(fc) 



^ (P(fc)+a^/A^.)'. (22) 



B. Combining multiple observations 

We now investigate the relative improvements in sen- 
sitivity from reobserving our sky patch, either using a 
different scan angle or the same scan angle as the first 
observation. We will refer to these strategies as "mini- 
mally cross-linked" and "non-cross-linked." 

As discussed in the Introduction, sub-orbital CMB 
missions can achieve a degree of cross-linking with 
azimuth-only scans by observing a patch of sky at dif- 
ferent times during the day. In particular, the minimally 
cross-linked strategy described above is achievable by ob- 
serving a sky field at the same elevation twice a day: once 
as the field rises, once as it sets. The maximum difference 
in scan angles available to a sub-orbital platform depends 
on the latitude of the observing platform: 



< 90° 



|lat|. 



(23) 



As noted earlier, no cross-linking is possible with 
azimuth-only scans from a Polar platform. 

If we reobserve our sky patch at the same scan angle 
(no cross-linking), the fraction of good modes will not 
change, but the noise variance in those modes will de- 
crease by a factor of 2. If we reobserve the patch at a 
different scan angle (minimal cross-linking) , we will have 
three classes of modes: 

1. Modes which satisfy equation [T51 in both observa- 
tions, 

2. Modes which satisfy equation [TH] in one observation 
but not the other, and 

3. Modes which satisfy equation [TH] in neither obser- 
vation. 
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FIG. 1: Graphical representation of the three classes of modes 
enumerated in section Hi Bl Left Panel: Phasor plot of modes 
for a single k — \k\. For two observations at scan angles 
9si and 9s2, modes which fall within 9c of either 9ai or 6s2 
but not both - indicated by the singly hatched regions - will 
have variance a^/Ns; modes which fall within 9c of both 9si 
and 9a2 - indicated by the doubly hatched region - will have 
variance a1,/2Ns\ modes which lie further than 9c from both 
Ssi and 9a2 (non-hatched regions) will have infinite variance. 
Right Panel: 2-d Fourier-space representation of these classes 
of modes at all k for ^si = 0° and 9^2 = 40°. The circle at 
the center has radius 2tt fiance/ Vs and indicates the modes that 
will have infinite variance no matter how many different scan 
angles are used. 



If we optimally combine the observations, these classes 
of modes will have noise variance a'^/2Ns, crJ/iVs, and 
oo, respectively (all modes will still have the same sam- 
ple variance P{k)). If the two scan angles are suffi- 
ciently different {\9si — Os2\ > 29c{k)), the two observa- 
tions will sample fully independent modes, and the frac- 
tion of modes falling into each noise variance category 
can be expressed as: 





iVi(fc) 


Ntotik) 


iVtot(fc) 




N2{k) 


iVtot(fc) 


iVtot(fc) 




Nsik) 


iVtot(fc) 


Ntotik) 





40c (fc) 



(24) 



1 



40c (fc) 



whereas if there is some overlap between modes sampled 
in the two observations (i.e., if j^si — 6*821 < 26'c(fc)), then 
the fraction of modes falling into each noise variance cat- 
egory is: 



iVi(fc) _ 20,{k) - %i - e,2\ 



Ntotik) 
N2ik) 

Ntotik) 
mjk) 

Ntotik) 



2% 



TT 
0s2| 



(25) 



= I - 



TT 

20c (fc) 



with the obvious limit that < Niik)/ Ntotik) < 1 in all 
cases and that J^i Niik) = Ntotik). Equations [Ml and[^ 
are represented graphically in figure [TJ 



To get the minimum-variance estimate of the CMB 
power spectrum from these three classes of modes, we 
estimate P(fc) using each class individually and make a 
weighted mean of these estimates using inverse- variance 
weighting based on equation 1151 In the noise-dominated 
regime (where a'^/Ns ^ Pik)), the variance on our mea- 
surement of the CMB power spectrum from these com- 
bined observations will be: 



[6Pik)y 



1 



[spih)r 

N^ik) , Niik) 



(26) 



.i</Ns)' 

= [SPnol/fik)]'- 



icrl/2N,)' 
N^ik) 



Niik) 



ANtotik) Ntotik) ^ 
where we have identified the prefactor in equation 



[6Pnol/fik)]' = 



^tot(fc) 



Pik) 



4iV2 A:/3kyAfc 



2A^, 

if al/N, » Pik) 



(27) 



as the variance our combined measurements would have 
if the instrument noise were equal to the high-frequency 
white-noise value at all temporal frequencies (i.e., with- 
out the 1// noise component). The ratio of variance in 
the presence of our toy- model 1// noise to the variance 
with no 1// noise (still assuming a1,/Ns ^ Pik)) evalu- 
ates to 



(5P(fc) 



1 2 



6Pnol/fik) 



(28) 



in the independent-modes case (where |0si— 0s2| > 
2e'c(fc)) and to 



6Pik) 



1 2 



JPnol/fik) 



26, 



7s2 



/2 



(29) 



< 



in the overlapping- modes case (where |0si — Os2 
20c(fc)). 

Equations [28] and [29] contain the potentially surpris- 
ing result that in the noise-dominated regime, the best 
constraints on P(fc) are obtained by reobserving the sky 
at the same scan angle (so that |0s2 — 0si| — 0). That 
is to say, the scan strategy without cross-linking out- 
performs the cross-linked strategy. It is trivial to show 
that in the sample-variance-dominated regime (where 
P(fc) ^ cr^/A's), one comes to the opposite conclusion, 
namely that the best constraints are obtained in the min- 
imally cross-linked strategy. In fact, the question being 
considered here is very similar to that of optimizing sky 
coverage vs. observing time per pixel in a CMB experi- 
ment. In our case the question is how to divide observing 
time between Fourier modes rather than spatial pixels, 
but the optimization calculation is the same. 
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C. Results for "real" 1// noise 

The approximation in equation [13] leads to a concise, 
illustrative result, but there is no reason we cannot do 
the calculation for noise with a true 1// spectrum (as 
defined in equation [T2| . For a single observation at scan 
angle 9s, our best inverse- variance weighted estimate of 
P{k) will have variance: 

(30) 



[sp{k)y 



2fc/skyAfc 
I.tt/2 

d{0k 




P{k) 



knee 



Optimally combining two observations at 9si and 9s2 
yields a power spectrum measurement with variance: 



[sp{k)Y 



2tt 



kfsky^k 




d9i, 



+ ^(1 + 91)11(1 +92) 



where 



ft 



2vr/kn 



kvs \cos{9k - 9s,i)\ 
and the parallel operator is defined such that 

i i "' ^ 

Xl\\X2 



1 

Xl 



1 

X2 



(32) 



(33) 



Once again, the relative performance of the non-cross- 
linked and minimally cross-linked strategies depends on 
the relative size of the noise and sample variance. Fig- 
ure shows the ratio of dP{k) to the value we would 
obtain with no 1// noise for the two scan strategies and 
two noise- vs. sample-variance regimes under considera- 
tion, using realistic values of Vg and /knoc- Though the 
difference in the two strategies is less dramatic than the 
factor of \/2 in the toy model case (see equations [28l and 
the behavior is qualitatively similar. In the noise- 
variance-dominated regime, the power spectrum sensitiv- 
ity is dominated by the best-measured modes, in which 
case it is advantageous to concentrate observing time on 
measuring a small number of modes very well; in the 
sample- variance-dominated regime, sensitivity is equal 
among all modes that are measured well enough to get 
below the sample variance, in which case it is better to 
spread the observing time around and measure as many 
modes as possible just well enough. 



D. Results for other choices of scan parameters 

While the most important conclusion to be drawn from 
figure [5] is the relative performance between the non- 
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FIG. 2: Ratio of 5P{k) to the no-1//- noise value for the non- 
cross-linked (solid line) and minimally cross-linked (dashed 
line) scan strategies in two regimes of noise variance vs. sam- 
(31) pie variance. Upper curves: P{k) = a^/Ns/iO. Lower curves: 
P{k) = 40 al/N,. In both cases, the observational param- 
_i eters are: /knee = O.lHz and Vs — l°/s. In the minimally 

cross-linked cases, ds2 ~ Oai + 90° . 



cross-linked and minimally cross-linked strategies, the 
amplitude of 5P{k) for both strategies compared to the 
case with no 1// noise is also of interest. Not surpris- 
ingly, in the noise-variance-dominated regime, the value 
of k at which the 1// noise begins to dominate 5P{k) is 
very close to what would naively expect, namely 



27r/i 



knee 



(34) 



For the particular values of /knee and Vs used in figure 
m (/knee — O.lHz and Vs — l°/s), which are realistic ap- 
proximations for the upcoming generation of CMB exper- 
iments, k{l/ f) = 36, indicating that for an experiment 
which meets these criteria, 1// noise will not limit obser- 
vations at least down to the £ ^ 100 range (recall that we 
are working in the flat-sky regime where £ ^ k). As for 
experiments which have different values of /knee and Vs, a 
quick glance at equations ED - ESI makes it clear that for 
the noise-variance-dominated case, the results in figure 
[2] can be scaled from the values for /knee and Vs used in 
the plot to arbitrary values by scaling the x-axis using 
equation [Ml 

The other key observing strategy design parameter for 
raster-scanned experiments is the size and geometry of 
the sky patch observed in one full raster. The results 
in figure [5] are unaffected by changes in sky patch size 
and geometry (which only come into equation [21] in the 
/sky term, which itself drops out when we normalize to 
the no-1// case), except in that the extent of the largest 
dimension of the patch puts a fundamental lower limit 
on the k- (or £-) modes that can be measured. 
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III. DISCUSSION 
A. Comparison with earlier work 

The results of equations [28] and [29] and figure [2] seem 
in direct contradiction to the conventional wisdom that 
cross-linking — however minimal — is vital to CMB 
power spectrum sensitivity in the presence of 1// noise. 
In particular, 0, T97] found an order-of-magnitude dif- 
ference in P"°'*'°(fc) for two scan strategies very much 
like our non-cross-linked and minimally cross-linked ones 
(cf. T97's "Serpentine" and "Fence" scans). However, 
the quantity plotted in T97 is effectively the noise power 
averaged over all modes of a given k (cf. T97, equation 
52), whereas the total instrument sensitivity to signal at 
a given k is proportional to the sum of the weight, or 
the inverse noise power, over all modes of that k. The 
key to reconciling our result and that of T97 is to under- 
stand that the excess noise power in the non-cross-linked 
or Serpentine scan is concentrated in a few easily iden- 
tifiable bad modes, and that the remaining good modes 
are actually measured better in this scan strategy than in 
the minimally cross-linked, or Fence scan. When the in- 
verse noise power is summed over all modes of a given k, 
it is the repeated measurement of the good modes that 
wins out over measuring new modes — at least in the 
noise- variance-dominated regime. 
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IV. CONCLUSION 

We have calculated the sensitivity to the CMB power 
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